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Abstract 

We describe a bijection between hexagonal lattice animals and a special type of square 
lattice animals. Using this bijection we adopt Maple packages that automatically gener- 
ates generating functions (and series expansions) for fixed square lattice lattice animals 
to that of fixed hexagonal animals on the two-dimensional hexagonal lattice confined to 
a strip < y < k, for arbitrary k. 
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1. Background 

Required reading: Symbol-Crunching with the Transfer-Matrix Method in Order to 
Count Skinny Physical Creatures [Zl]. 

In [Zl], Zeilberger used finite transfer matrix method and developed two Maple pack- 
ages ANIMALS and FreeANIMALS to count square lattice animals confined to a strip. Here 
we define a bijection between hexagonal lattice animals and a special class of square lat- 
tice animals and adopt the packages HexANIMALS and HexaFreeANIMALS accompanying 
[Zl] to enumerate hexagonal lattice animals. 

Definition 1.1 A hexagonal lattice animal (hexagonal polyominoes hence forth) on 
the hexagonal lattice is an edge-connected set of lattice cells on the hexagonal lattice. 
Two animals are equivalent if they are translations of each other. 



1 This article is accompanied by two Maple packages, HexANIMALS and HexaFreeANIMALS, that can 
be downloaded from http : //www . people . vcu . edu/ ~mapagodu/ 



Figure 1 shows the three hexagonal polyominoes comprised of two hexagons. 




Figure 1: The three hexagonal polyominoes with two hexagons. 



The number of nonequivalent hexagonal polyominoes with n cells, a(n), is given by 
Sloane's sequence A001207 for n < 35 [NJAS]. For example, a(l) = 1, a(2) = 3, 
a(3) = 11, a(4) = 44, a(5) = 186, a(6) = 814, a(7) = 3652, a(8) = 16689, a(9) = 77,359, 
and a(10) = 362,671. 

A hexagonal polyominoes in which each column contains at most k contiguous blocks 
of cells is referred to as a /c-board polyominoes. Figure 2 shows examples of one and 
two-board polyominoes with 7 hexagons. The number of nonequivalent one-board poly- 
ominoes with n hexagonal cells, b(n), is given by Sloane's sequence A059716 for n < 24 
[GV]. For example, 6(1) = 1, 6(2) = 3, 6(3) = 11, 6(4) = 42, 6(5) = 162, 6(6) = 626, 
6(7) = 2419, 6(8) = 9346, 6(9) = 36, 106, and 6(10) = 139,483. 



Figure 2: A one-board hexagonal polyominoes (a) and a two-board hexagonal polyomi- 
noes (b) with 7 hexagons. 

The generating function that enumerates one-board hexagonal polyominoes is com- 
puted in [K], which can also be computed using the Maple package LEGO in [Z0] by 
taking p(a, 6) = a + 6 (for definition of p(a, 6) refer to [Z0]). To the best of the authors' 
knowledge, &-board hexagonal polyominoes have not been enumerated for k > 2. In 
this article, in addition to computing generating functions (and series expansions) for 
hexagonal lattice animals that fit into a prescribed but arbitrary height, we also compute 
the first 12 terms of the sequence that enumerate the number of board-pair-hexagonal 
polyominoes, the analog of square board-pair-polyominoes computed in [Z2]. Fisrt we 
recall some definitions and introduce the essential details of the Maple implementation 
of the transfer-matrix method [S] from [Zl]. 





(a) 



(b) 
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2. The Transfer-Matrix Method 



Let G(V, E) be a vertex weighted directed graph consisting a finite set of vertices, V, and 
a finite set of directed edges, E. A path P in G is a sequence v±,ei, e 2 , ■ ■ ■ , e m , v m 
where Vi G V(l < i <m) and e« G £"(1 < i < m) is an edge from to Vi. The weight 
of a path P is the sum of the weights of all vertices participating in the path P. In this 
case Wt(P) = Wt(vi) + Wt(v 2 ) + ... + Wt(v m ). 

For T a subset of V, we want to compute the weight enumerator (generating function), 

oo 

F(z)=5>,-*»" 

where Oj is the number of paths with weight j that starts at any vertex in V and ends at 
a vertex in T. For ueV, let F v (z) be the generating function that enumerates all paths 
starting at vertex v and ends at a vertex in V. Then, clearly, 

F v (z)= Yl tWt(P) C 1 ) 

paths P 
mit(P)=t), /m(P)eT 

where init(P) is the initial vertex of P and fin(P) is the terminal vertex of P. Let 
iV(i>) be the set of vertices adjacent to v. Then, {F v (z)\v G V} satisfies the following 
\V\ (where \V\ is the number of vertices) system of equations in |V| unknowns, namely 
{F v (z)\veV}, 

F v (z) = l T (v)z wt W + z wt W n(v,u)F v (z) (2) 

ueN(v) 

where n(v, u) is the number of paths from v to u and 

l( z )-f li f zeT (o\ 

T 1^0 otherwise 

Since each F u (z) exists as a power series in z, our system has a solution. Once we 
solve for this sytem, the answer to our original problem is then 

F(z) = J2F v (z) 
i>eT 
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As an example consider the following vertex-weighted directed graph on four vertices 
with weights as shown next to the vertices. 



w,2 




u,l v,2 z,l 



Figure 3: Example of vertex weighted directed graph on four vertices 

In this example, the generating function that enumerates all paths according to their 
weight is found by first finding the four generating functions {F u (z), F v (z), F w (z), F t (z)} 
for all paths starting at vertex u, v, w, and t respectively. These four generating functions 
are related by 

ft(z) = z, f w {z) = f v (z) = zf u (z) = . 

1 — z° 

Solving these system and adding together we get the required generating function that 
enumerates all paths according to their weight as 

„, , 2z + 2z 2 + z 3 + 3z 4 + 3z 5 + z 6 
F(z) = . 

For the hexagonal polyominoes case, we follow the structure of [Zl] and automate the 
problem and use Maple do the hard part. First, we recall the following definitions from 
[Zl]: 

Definition 2.1 A Combinatorial Markov Process is a six tuple (V, E, init, fin, Start, 
Finish), where V is a finite set of vertices, E is the set of directed edges, and init, fin : 
E i — ^ V are functions that assign the initial and terminal vertex to an edge, respectively. 

Definition 2.2 A vertex Weighted Combinatorial Markov Process is a seven tuple 
(V, E, init, fin, Start, Finish, wt), where (V, E, init, fin, Start, Finish) is a Combinatorial 
Markov Process and wt : V(E) i— > Z + is a function that assigns weights to vertices. 

In Maple, we represent a Combinatorial Markov Process with vertices V — {1,2, ... ,n} 
and no multiple edges as a four tuple [S,T,ListOf OutGoingNeighbors ,ListOf Weights] 
where S,T C V, and ListOf OutGoingNeighbors and ListOf Weights are lists of length 
n whose i th element is the set of vertices adjacent to vertex % and the weight of vertex i, 
respectively. 
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Similarly, the Combinatorial Markov Process for a multiple edge graph is a four tuple 
[S,T,ListOfOutGoingNeibors,ListOf Weights] where S, T, and ListOf Weights are as 
previously defined and ListOf OutGoingNeighbors is a list of length n whose i th element 
is a multi-set j^j™ 2 • • • fk k represented in the form {[ji, mi], \j 2 , m 2 ], . . . , \jk, mk]}, which 
means that from vertex % there are m\ edges to vertex ji, m 2 edges to vertex j 2 , and so 
on. 

If we successfully model hexagonal polyominoes as a weighted Combinatorial Markov 
Process as it is done for the squre lattice animals in [Zl] then we can employ the Maple 
package MARKOV given in [Zl] to automatically find the generating function and series 
expansion of the generating function. The next section describes how weighted Combi- 
natorial Process for hexagonal polyominoes is created. 

3. Maple Representation of hexagonal polyominoes 

We could develop a grammar that describes height-restricted hexanimals from first prin- 
ciples and apply the same methods as Zeilberger [Zl] to encode the grammar as a Com- 
binatorial Markov Process; however, this would result in rewriting much of the code that 
was already created for the square lattice animals. Our approach here is to define a 
bijection between hexagonal lattice animals and a special class of square lattice animals 
and adopt the Maple codes ANIMALS and FreeANIMALS in [Zl]. In the remaining 
sections, we describe this processes. 

We overlay a hexagonal polyominoe on the square lattice so that its leftmost vertex 
(vertices of the square contained inside the hexagonal) lie on the line x = and its bottom 
most cell lie on the line y = as shown in Fig.4(a). We then construct a polyominoes 
comprised of the square cells which have diagonally opposing corners falling on the same 
hexagonal cell as shown in Fig. 4(b). 

As described in detail by Zeilberger [Zl], we represent a polyomino by the set of 
coordinates of the bottomleft corner of each of its cells. We further encode a polyomino 
as a word in the alphabet consisting of the non-empty subsets of {0, 1, k} where k is 
the maximum y coordinate in the set representation of the polyomino. For example, the 
animal in Fig. 4 (b) is represented by the set 

{(0,2), (0,3), (1,1), (1,2), (1,5), (1,6), (2,0), (2,1), (2,2), (2, 3), (2,4), (2, 5), (3,1), (3,2)} 
and encoded as the word (in interval notation) 

{[2, 3]}, {[1,2], [5, 6]}, {[0,5]}, {[1,2]}. 

There are two important properties of the polyomino words that encode hexanimals. 
First, every interval has even size (i.e., it represents a contiguous block of an even number 
of cells); this is obvious from the mapping which associates each hexagonal cell with a 
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Figure 4: Mapping between a hexagonal lattice animals and a square lattice animals 



pair of square cells. Second, within a letter, the starting position of every interval has the 
same parity (odd or even), and adjacent letters have different parity; this is again obvious 
from the mapping of each hexagonal cell to a pair of square cells and the vertical shift of 
adjacent hexagonal columns. Note, however, that the parity of the leftmost column may 
be either odd or even depending on the hexanimal. 

We will call the class of square polyominoes to which hexagonal polyominoes are 
mapped parity polyominoes. It is clear that the mapping between hexagonal animals and 
parity polyominoes is unique and reversible, and thus defines a bijection. As a result 
of this bijection, we can count hexagonal animals by counting parity polyominoes. Zeil- 
berger's code counts general polyominoes (a much larger class than parity polyominoes), 
so in order to use it to count parity polyominoes, we need to restrict the alphabet. In 
the following sections of this paper, we describe how to modify Zeilberger's code to count 
parity polyominoes while minimizing the amount of changes required; we assume that 
the reader is familiar with [Zl]. 



4. Counting Globally Skinny Hexagonal Animals 



Zeilberger's ANIMALS package uses the transfer-matrix method [S] to count polyominoes 
whose set representation has y coordinates in the range < y < k for some arbitrary k. 
Before the transfer-matrix method can be applied, a grammar needs to be constructed 
to describe how to combine letters into words that encode polyominoes. A key routine 
in the generation of this grammar is PreLef tLetters (a,b) which returns the set of 
all possible subsets of the integers in the range [a, b] written in interval notation. For 
example, PreLef tLetters (0,2) returns the following subsets: 

{{}, {[0,0]}, {[0,1]}, {[1,1]}, {[0,2]}, {[1,2]}, {[2,2]}, {[0,0], [2,2]}} 
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For parity polyominoes, we are only interested in those subsets whose intervals have 
even size and the same parity for the starting position. Rather than rewriting 
PreLef tLetters from scratch, we can use a trick in Maple to subclass PreLef tLetters 
and alter its parameters without having to modify the original PreLef tLetters code: 

read ANIMALS: 

origPreLef tLetters := subs (PreLef tLetters=origPreLef tLetters, 

eval (PreLef tLetters) ) : 

PreLef tLetters := proc(a, b) 

local halfS, S, i, PreLet, PreLetl: 

halfS := origPreLeftLetters(0, floor((b-a+l)/2)-l) minus {{}}: 
S := {{}}: 

for i from 1 to nops(halfS) do 

PreLet : = {seq( [a+2*half S [i] [j] [1] , a+2*half S [i] [j] [2] +1] , 

j=l. .nops(halfS[i]))}: 
PreLetl := {seq( [PreLet [j] [1] +1 , PreLet [j] [2] +1] , 

j=l. .nops (PreLet))}: 
S := S union {PreLet}: 

if max(seq(PreLetl [j] [2] , j=l . .nops (PreLetl) ) ) <= b then 
S := S union -[PreLetl}: 

fi: 

od: 
S: 

end: 

The first line loads Zeilberger's original ANIMALS package. The second line makes a 
copy of PreLef tLetters; the subs function is designed to modify PreLef tLetters so 
that any recursive calls are made to the copy rather than the original. The remaining 
lines redefine PreLef tLetters so that it produces the type of intervals that are required 
by parity polyominoes. The new PreLef tLetters uses the original PreLef tLetters 
to generate a set of subsets (in interval notation) for half the required range (stored in 
halfS). The for loop goes through each interval and doubles its size; this results in 
intervals that have even size and even parity (stored in PreLet). To create intervals with 
even size and odd parity, each interval in PreLet is shifted down by one, and the result 
is stored in PreLet2. The final condition checks that the shifted intervals do not go out 
of the range [a, b]. 

For example, PreLef tLetters (0,5) would call origPreLef tLetters (0,2) and get 

the subsets shown above. The for loop would then modify these subsets to generate the 
following result which represents the parity polyomino preletters in the range [0,5]: 
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{{},{[0, !]},{[!, 2]}, {[0,3]}, {[1,4]}, {[2, 3]}, {[3, 4]}, {[0,5]}, {[2, 5]}, {[4, 5]}, {[0,1], [4, 5]}} 



Once the preletters are generated, the grammar is extended by determining which of 
the preletters are valid extensions of existing polyomino words; the ANIMALS routine that 
performs this check is PreLetToLet. For parity polyominoes, not only must the preletter 
satisfy the same connectivity constraints as for general polyominoes, but the parity of 
the preletter's interval starting positions must be different than the letter it follows. As 
we did with PreLef tLetters, a simple extension to the original code is that is required: 

origPreLetToLet := subs(PreLetToLet=origPreLetToLet , 

eval (PreLetToLet) ) : 

PreLetToLet := proc(Let, PreLet) 

if PreLet <> O and Let[l] [1] [1] mod 2 = PreLet [1] [1] mod 2 then 
0: 

else 

origPreLetToLet (Let , PreLet): 

fi: 

end: 

Since we know that every interval in a parity polyomino letter has the same starting 
position parity, it suffices to check whether or not the starting positions of the first 
intervals have differing parities; if that is the case, PreLetToLet calls the original ANIMALS 
routine to ensure that the connectivity constraints are satisfied. 

In order to weight the transfer-matrix vertices so that polyominoes with a specific 
number of cells can be enumerated, the ANIMALS routines call Weight (Let) to determine 
the number of cells in the letter Let. For enumerating hexanimals, we want to weight 
the transfer-matrix vertices with the number of hexagonal cells. Since each hexagonal 
cell is represented by two polyomino cells, a letter's weight in hexagonal cells is half its 
weight in polyomino cells: 

origWeight := subs(Weight=origWeight , 

eval (Weight) ) : 

Weight := proc(Let) 

floor (origWeight (Let) /2) : 

end: 
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One final change to ANIMALS is needed in order to enumerate hexanimals; the routine 
Khaya(L) that computes the number of polyominoes with < L cells relies on symmetry 
that does not exist in the hexagonal lattice. Since a hexanimal with L cells must map to 
a parity polyomino whose set representation has a maximum y coordinate of 2L, we can 
implement Khaya(L) by a call to GFseries: 

Khaya := proc(L) 

[GFseries(2*L, L)] : 

end: 



5. A User's Manual for HexANIMALS 

The modifications to ANIMALS that are described in the previous section are contained 
in the package HexANIMALS. To run HexANIMALS, download it and ANIMALS from the 
first authors' web site to a local directory and start Maple. Once in Maple, type: read 
HexANIMALS and follow the on-line help. 

Excluding Khaya, all the generating function and series expansion routines in ANIMALS 
remain unchanged. GF(n,s) computes the generating function for hexanimals embed- 
ded in the square lattice of height n (i.e., the height of the square lattice is comprised 
of n cells). GFseries (n,L) computes the list of length L whose k th term is the num- 
ber of hexanimals with k cells embedded in the square lattice of height n. Khaya (L) 
computes the list of length L whose k th term is the number of hexanimals with k cells. 
Gf and Gf series are the analogs of GF and GFseries for hexanimals whose polyomino 
embeddings have height exactly n. 

As an example of HexANIMALS usage, the call GF(l,s) returns since no hexanimal 
can be embedded in a square lattice of height 1. GF(2,s) returns s since there is exactly 
one hexanimal (the single cell) that can be embedded in a square lattice of height 2. 
GFseries (3, 5) returns [1,2,2,2,2] since the hexanimals in the square lattice of height 
3 with more than one hexagonal cell are chains, and each chain has two orientations 
depending on the parity of the starting cell. In Table [1] we compute the generating 
function for hexanimals embedded in a square lattice of height n < 7 along with the first 
few terms in the expansion by hexagonal cell count. Once n > 13, it takes too long to 
compute GF(n,s) exactly, but one can go much further with GFseries(n,L). 

6. Counting Locally Skinny Hexagonal Animals 

If we relax the condition that the entire hexagonal polyominoes must fit within a square 
lattice of fixed height, and only require that each column have height less than some ar- 
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n 


F(z) 




1 





0,0,0,... 


2 




1,0,0,... 


3 


2(1 + 2) 
1-2 


1,2,2,... 


4 


2(1+2) 


1,3,6,11,19,32,53,87 




5 


2(2^+22^+2 + 1 


1,3,10,25,61,142,323,723 


( z 3 + 2z 2 +z-l)(z 2 +z-l) 


6 


2(l + 2 + 22 2 +72 b +42' :1 +82 4 -32 b -42 Y +32 8 +2 1 0+42 M ) 


1,3,11,37,111,320,896 


(32 6 +22 7 -l +22 2 +2+52 3 )(2 3 +22 2 +2-l) 



Table 1: Generating functions for hexanimals polyominoes embedded in the square lattice 
of height n and their series expansion by hexagonal cell count. 

bitrary constant, then we get what Zeilberger calls locally skinny hexagonal polyominoes 
(as opposed to globally skinny hexagonal polyominoes). 

The code for computing the generating function and series expansions for locally 
skinny polyominoes is contained in Zeilberger's FreeANIMALS package. Since the changes 
to FreeANIMALS that are required in order to enumerate hexagonal lattice animals are 
virtually the same as those made to ANIMALS, we omit the details that were discussed in 
the previous sections and only discuss the major differences. 

FreeANIMALS creates a Combinatorial Markov Process with multiple edges between 
letters. Each edge represents one of the possible vertical offsets of the letters. In contrast, 
ANIMALS creates a Combinatorial Markov Process with only one edge between letters; 
thus, each letter represents a particular arrangement of cells with an intrinsic offset. As 
a result, PreLef tLetters in HexANIMALS had to create the two possible parity offsets for 
each letter whereas PreLef tLetters in HexaFreeANIMALS needs to only create a single 
instance of the letter and the code that generates the Combinatorial Markov Process will 
create its own instances of the shifted letter. The resulting PreLef tLetters routine is 
similar to the PreLef tLetters routine in HexANIMALS except that the parity shifting 
code is omitted: 

read FreeANIMALS: 

origPreLef tLetters := subs (PreLef tLetters=origPreLef tLetters , 

eval (PreLef tLetters) ) : 

PreLef tLetters := proc(a, b) 
local halfS, S, i, PreLet: 

halfS := origPreLeftLetters(0, floor((b-a+l)/2)-l) minus {{}}: 
S := {{}}: 

for i from 1 to nops(halfS) do 
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PreLet := {seq( [a+2*half S [i] [j] [1] , 

a+2*halfS[i] [j] [2]+l] , j=l. .nops (half S [i] ) )} : 
S := S union {PreLet}: 

od: 
S: 

end: 

In addition, FreeANIMALS contains the routine PreLef tLettersk(a,b,k) which com- 
putes all possible subsets of the integers in the range [a, b] that are represented by exactly 
k intervals. PreLef tLettersk is used to compute &-board polyominoes. For example, 
PreLef tLettersk (0,4, 2) returns the following subsets: 

{{[0,0], [2, 2]}, {[0,0], [2, 3]}, {[0,0], [3, 3]}, {[0,1], [3, 3]}} 

The modifications to convert PreLef tLettersk for use with HexaFreeANIMALS are 
analogous to the changes made to PreLef tLetters as shown in the following code: 

PreLef tLettersk := proc(a, b, k) 
local halfS, S, i, PreLet: 

halfS := origPreLeftLettersk(0, f loor((b-a+l)/2)-l, k) minus {{}}: 
S := {{}}: 

for i from 1 to nops (half S) do 

PreLet := {seq( [a+2*half S [i] [j] [1] , 

a+2*halfS[i] [j] [2]+l] , j=l. .nops (half S[i] ))}: 
S := S union {PreLet}: 

od: 
S: 

end: 

The only other significant difference between HexANIMALS and HexFreeANIMALS is in 
the specification of height restrictions. When considering globally skinny hexanimals, the 
height restrictions were defined in terms of the square lattice in which the hexanimals 
were embedded; this was because the shifting of adjacent columns in hexanimals makes 
the definition of the height of a hexanimal in terms of hexagonal cells difficult. In the 
case of locally skinny hexanimals, the height restrictions apply to individual columns 
and are well-defined in terms of hexagonal cells. Since the underlying algorithms of 
HexaFreeANIMALS operate on parity polyominoes, the height restrictions for hexanimal 
columns are simply doubled when applied to parity polyominoes as shown in the following 
code: 
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origgf := subs (gf =origgf , eval(gf)): 

gf := proc(n, s) 

origgf (2*n, s) : 

end: 

origgf Series := subs(gfSeries=origgf Series, eval (gf Series) ) : 

gfSeries := proc(n, L) 

origgf Series (2*n, L): 

end: 

origgf List := subs(gf List=origgf List , eval (gf List) ) : 

gfList := proc(resh, s) 

origgf List ( [seq(2*resh[i] , i=l . .nops (resh) )] , s) : 

end: 

origgf SeriesList := subs (gf SeriesList=origgf SeriesList , 

eval (gf SeriesList) ) : 

gf SeriesList := proc(resh, L) 

origgf SeriesList ( [seq(2*resh [i] , i=l . .nops (resh) )] , L) : 

end: 



7. A User's Manual for the HexaFree ANIMALS 



The modifications to FreeANIMALS are contained in the package HexaFreeANIMALS. 
gf (n,s) computes the generating function for hexagonal polyominoes whose columns 
span < n hexagonal cells, gf Series (n,L) computes the list of length L whose k th term 
is the number of hexagonal polyominoes with k cells whose columns span < n hexagonal 
cells. gfList (List , s) computes the generating function for polyominoes whose /c-board 
columns span < List[k] hexagonal cells, gf SeriesList (List ,L) computes the series ex- 
pansion of 
gfList (List , s) up to L terms. 

As an example of the usage of HexaFreeANIMALS, the call gfList ([7, 5] ,s) would 
compute the generating function for polyominoes whose 1-board columns span < 7 hexag- 
onal cells and whose 2-board columns span < 5 hexagonal cells, gf SeriesList ( [24] , 24) 
computes Sloane's sequence A0059716 [NJAS] of 1-board polyominoes. 

The call gf SeriesList ( [12 , 12] ,12) computes the first 12 terms of the previously 
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unknown sequence 



1, 3, 11, 44, 186, 814, 3648, 16611, 76437, 354112, 1647344, 7682237 

that enumerates 2-board hexanimals with a given hexagonal cell count, the analog of 
Sloane's sequence A001170 for board-pair-pile polyominoes. 



8. Conclusion 



We have described the Maple packages HexANIMALS and HexaFreeANIMALS. HexANIMALS 
enumerates globally skinny hexagonal polyominoes in which each polyominoes is em- 
bedded in a horizontal strip of prescribed height. HexaFreeANIMALS enumerates locally 
skinny polyominoes in which the total height of each hexagonal polyominoes is unbounded 
but whose columns are composed of a bounded number of hexagonal cells. In addition, we 
also demonstrated used HexaFreeANIMALS to compute the first few terms of the number 
of board-pair-hexagonal animals. 
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